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Abstract 

This work presents concisely the results obtained from the analysis of the two-point function of the gauge 
field in the SU(2 ) and SU{2) x U{ 1) gauge theories, in the Landau gauge, coupled to a scalar Higgs field in 
the fundamental or adjoint representation. Non-perturbative effects are considered by taking into account the 
Gribov ambiguity. In general, in both Yang-Mills models the gluon propagator has non-trivial contributions of 
physical and non-physical modes, which clearly depends on the group representation of the Higgs field. These 
results were presented during the 4 th Winter Workshop on Non-perturbative Quantum Field Theory, which 
took place in Sophia-Antipolis - France. 


1 Introduction 

The understanding of non-perturbative aspects of non-Abelian gauge theories is one of the main challenging 
problems in quantum field theories. As an example, we may quote the transition between the confined and the 
Higgs regime in a Yang-Mills theory coupled to a scalar Higgs fields. See refs.[l, 2, 3] for analytical investigations 
and [4, 5, 6, 7, 8, 9] for results obtained through numerical lattice simulations. 

Non-perturbative effects can be accounted perturbatively by considering ambiguities in the gauge fixing 
process, first noted by Gribov in [10]. These ambiguities, also referred to as Gribov copies, are unavoidably present 
in the Landau gauge, since the (Hermitian) Faddeev-Popov operator admits the existence of zero modes. A brief 
introduction to the Gribov issue is presented in section 2, based on the works [10, 11, 12]. With that being said, 
the analysis of Yang-Mills theories coupled to a scalar Higgs field, in its fundamental and adjoint representation, 
while taking into account non-perturbative effects arising from the Gribov issue, becomes of great aid in the 
comprehension of the confined/Higgs transition effect. More specifically, the SU(2) and SU(2) x U(l) Yang- 
Mills models are analysed in this work, in the 3- and 4-dimensional cases. Some general features and definitions 
concerning the Yang-Mills + Higgs field are worked out in section 3. 

Part of the dynamics of the system can be read off from the structure of the gluon propagator, which 
is a function of the parameter space of each model and reflects the existence of Gribov’s issue. The Higgs and 
confinement regimes are detected as follows: when the gauge field displays a Yukawa type propagator the system 
is said to be in the Higgs regime ; for a propagator of the Gribov type the system is in the confined regime. The 
confinement interpretation follows from the Gribov propagator that lacks the usual physical particle interpretation, 
since it exhibits complex conjugate poles. As such, gluons do not belong to the physical spectrum of the theory, 
which is compatible with the effect of confinement [13, 14]. Note that the “phase structure” in each theory turns 
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out to be dependent on the group representation of the Higgs field. The analysis of each model, SU( 2) and 
SXJ(2) x 17(1), are carried out in sections 4 and 5, respectively. In the last section 6 our conclusions are presented. 
The content of the present work is a summary of the results found by the authors in the works [15,16, 17, 18]. 


2 A brief introduction to Gribov’s issue 

This section is organized in a way to provide a brief introduction to Gribov’s ambiguities by following his 
seminal work [10]. We do not intend to provide the final word on this matter and as such we would refer to [12] 
for a more complete and pedagogical reference. The Faddeev-Popov quantization procedure is reproduced (not in 
full detail) in the subsection 2.1, while Gribov’s problem will be introduced and implemented on the path integral 
in the two subsequent subsections, 2.2 and 2.3. In the last subsection 2.4, we should present one of the most 
important outcomes of the Gribov’s issue: a possible interpretation of gluon confinement, which is encoded in the 
poles of the gluon propagators. We should say, to clarify matters, that our present work concerns computations up 
to one-loop in perturbation theory. 

2.1 Quantization of non-Abelian gauge field 

The gauge invariant action of a non-Abelian gauge field, or the Yang-Mills (YM) action, is given by 

Sym = J d d x X F; v F; v , (i) 

with F® v = + gf abc A b ^Al being the field strength tensor. The action (1) enjoys the feature of being 

symmetric under gauge transformation, which is defined for the gauge field as 

a ; = U'ApU-'-UHdrU), ( 2 ) 

with U{x) £ SU(N) and N being the number of colours. A geometrical representation of this symmetry can be 
seen in Figure 1, where each orbit, representing equivalent gauge fields, is crossed by a gauge curve F. It means 
that, the YM action is invariant under transformations that keep the transformed gauge field on the same gauge 
orbit. 

As is well known, that gauge invariance of the action induces inconsistencies in the quantization of the 
gauge field reflecting the existence of infinite physically equivalent configurations. To get rid of those spurious 
configurations from the system one has to fix the gauge, which is, in the geometrical view, to choose a convenient 
curve F that crosses only once each gauge orbit. In the path integral formalism, the gauge fixing procedure is 
carried out by the Faddeev-Popov procedure. 



Figure 1: Gauge orbits of a system with rotational symmetry in a plane and a function F which picks one representative 
from each gauge orbit. 
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The (Euclidean) gauge fixed partition function reads 


Z YM {J) = J [(LA] e-^+f^W , 

where Jj. denotes the path integral restricted to the curve F, which can be recast in the form 


( 3 ) 


Z YM (J) = V J[dA] Ar(A)6(F(A))e- s ™+f dxja ^ . (4) 

The V factor accounts for the (infinite) orbit’s volume, while Ajr(A) stands for the Jacobian of infinitesimal gauge 
transformations 

4 = A„ - D,e(x) , (5) 

which is needed since we are working with the gauge transformed integration measure. The 0(x) = 9 a (x)r a stands 

for the infinitesimal gauge transformation parameter, while r a are the SU(N) generators. <5(.F(A)) is the delta 
function ensuring the gauge condition F(A) = /. We should emphasize that the Jacobian of a given transformation 
(as the infinitesimal gauge transformation) is defined as the absolute value of the determinant of the derivative — 
with respect to the transformation parameter ( 6 ) — of the transformed field. In our gauge fixing case we have 

, ( 6 ) 

J r {A ’)=0 

with 


Ajf(A) = | det M a b(x,y)\ with M ab (x,y) = 


dF^A'^x)) 

50 b (y) 


F a {^{x)) = d„A'“(x) - , (7) 

being the gauge condition and A' given by (5). b a stands for the Lagrange multiplier: an arbitrary scalar field 
(the Nakanishi-Lautrup field). With this we have the Faddeev-Popov operator A4 ab (x,y), which reads 

Mab{x,y ) = -d^Df5{y-x)\ HA)=Q . ( 8 ) 

Following the standard procedure by Faddeev-Popov to implement the gauge fixing, one ends up with 


Z(J) = / [dA] [dc] [dc] exp 


—S- 


YM 


d.T 


c-d,D; b c b 


- 




Sgf 


(9) 


where S g f stands for the gauge fixing term. The ghost fields (c, c) are anti-commuting fields with ghost number 
+1 and —1, respectively. At the end we have the Faddeev-Popov action given by iSfp = Sym + Sgf. 

Let us emphasize two important assumptions made in the process to obtain (9): 

• The gauge condition F a is said to pick up only one field configuration from each gauge orbit, representing 
the physical equivalent configurations related by gauge transformations; 

• The determinant of A4 ab (x,y ), present in the Jacobian definition (6), is supposed to be always positive. 

Gribov showed in [10] that those assumptions are not always true. Furthermore, he showed that the failure of the 
first statement is closely related to the failure of the second one and, besides that, he pointed out the existence of 
zero mode configurations related to the FP operator (8). Those subtleties give rise to the Gribov issue, summarized 
in the next subsection. 


2.2 Gribov’s issue 

As stated before, one of the most important hypotheses is that once the gauge fixing condition is chosen, 
F = 0, for a given gauge field configuration A' , one should be able to find out only one gauge field configuration A M 
related to A' through a gauge transformation that fulfils the gauge condition. This situation is represented in Figure 
2 by the curve L. It means that we cannot find gauge-equivalent configurations, related by gauge transformation, 
satisfying the gauge condition (7), which is represented by the curve L'. Besides that, we cannot find a situation 
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Figure 2: The gauge condition curve can cross each orbit of equivalence once, more then once and at no point at all. The 
horizontal axes denotes the transversal gluon propagator, while the vertical axis represents the longitudinal one 

[ 0 ]. 


where no one field configuration on a gauge orbit satisfy the gauge condition (7), which is represented by the curve 

L". 

While there is no known example where the curve L" does exist, the situation described by curve the L' 
is quite typical in non-Abelian gauge theories and, therefore, it is worth analysing [10, 11, 12]. To this end, let us 
consider two close gauge-equivalent configurations: A' and A M related by an infinitesimal gauge transformation 
(5) and satisfying the Landau gauge condition (covered for £ —» 0). That is, 

A' m = UA„Ut- , d = 0 & d M A' M = 0. (10) 

Then, we are led to 


-w* = °> (H) 

with the following covariant derivative Z2“ fc = 9 M <5 a6 — gf abc A c Equation (11) does point to the existence of zero 
eigenvalues of the FP operator. The conclusion is that if there are (at least) two close equivalent fields satisfying 
the Landau gauge, defining gauge copies, then the Faddeev-Popov operator has zero modes (eigenstates associated 
to zero eigenvalues). Note that for the Abelian case the equation reduces to the Laplace equation, 

d 2 a = 0 . (12) 

Since it defines plane waves, thus not normalisable, we are not considering this case, restricting ourselves to the 
analysis of fields that vanish at infinity. It is clear one needs to study the FP operator’s eigenvalue. To that end, 
suppose the follow eigenvalue equation 


- d^D^ip = eip . 

(13) 

Clearly, if is small enough then we have 

- dfy = ^ , 

(14) 


which is solvable only for positive e, since in the momentum space we have p 2 = e > 0. Then, for small enough 
field configurations we do not have zero mode issues anymore. Otherwise, if turns out to be large, but not 
too large, we reach the solution e = 0, showing the existence of zero modes. For even large amplitudes of A M the 
eigenvalue turns to be negative; if it keeps growing then e reaches zero again, and so on. 

In his work [10] Gribov showed that the domain of functional integration should be restricted to the first 
region, named “first Gribov region” 52, where e > 0, in order to avoid gauge copies. It is known that this region is 
actually not free of copies. It has being analytically motivated that the existence of those copies inside 52 do not 
influence physical results (see [19] and references therein). Furthermore, up to now there is no way to implement 
analytically the restriction of the partition function to the region free of copies (known as Fundamental Modular 
Region, A). 

The first Gribov region 52 enjoys some useful properties that have mathematical proof only in the Landau 
gauge, which justifies our interest in this gauge (see [12] and references therein). Namely, 

• For every field configuration infinitesimally close to the border <552 and belonging to the region immediately 
out side 52 (called ffy), there exist a gauge-equivalent configuration belonging to 52 and infinitesimally close 
to the border <552 as well [10]. It was also proven that every gauge orbit intersects the first Gribov region 52 
[ 20 , 21 ], 
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• The Gribov region is convex [22]. This means that for two gluon fields A 1 and A ^ belonging to the Gribov 
region, also the gluon field A M = a A* + /3 A 2 with a, /3 > 0 and a + /3 = 1, is inside the Gribov region. 

• One may also show that the Gribov region is bounded in every direction [22]. 

For more details concerning (proofs of) properties of Gribov regions see [12] and references therein. 

2.3 Implementing the restriction to the first Gribov region 0 

As was discussed before, our aim is to restrict the domain of the functional integration to the region where 
the FP operator has positive eigenvalues. Therefore, let us define the first Gribov region as the region where the 
FP operator is positive definite. Namely, 

n = {A;, d^A; = 0, M ab > 0} , (15) 


where A4 ab is the FP operator, 

M ab (x, y) = - d,D; b 6(x - y) = (-d 2 ^ + d,f abc A^) 6(x - y) , 
and the condition on this operator defines a positive definite operator, given by 


dxdyu a (x)M ab (x,y)u; b (y) 


> 0 


(16) 


(17) 


Since the FP operator is the inverse of the ghost propagator, the ghost propagator plays a central role. With this 
being said, we are going to compute the ghost propagator till the first loop order, by following [10]. The effective 
implementation of the restriction of the functional integration will be done by imposing the no-pole condition, 
which will become clear in what follows. The computation of the ghost propagator can be done by making use of 
the partition function (9). Let us denote the ghost propagator as 

(c a (p)c b (-p)) = G(k 2 1 A) ab . (18) 


Treating the gauge field as a classical external field we get 


G(k 2 ,A) 


whereby 


jy 2 j fiabG{k i-'tyab 

p (1 + <r(k,A)) , 


1,11 N g 2 

¥ + v¥n 2 -1 


d d q 

(27r) fl 


Ai(-q)Ai(q) 


(k - q)y.g v 

(k - q) 2 


er(fc, A) 


1 1 Ng 2 
Vk 2 N 2 - 1 


d d q 

(2ir) d 


Al(-q)Ai(q) 


(k ~ q)y.k v 
(■k - q ) 2 


The ghost propagator (19) can be perturbatively approximated by 


(19) 


( 20 ) 


G(k 2 ,A) 


1 1 
k 2 (1 — a(k, A)) 


( 21 ) 


whereby we can see that for a(k , A) < 1 the domain of integration is safely restricted to f2, characterising the 
no-pole condition. It can be proved that a(k,A) is a decreasing function of k , which means that the largest value 
of a is obtained at k = 0. Thus, if the condition a( 0, A) < 1 is ensured then it is completely safe for any non zero 
value of k. Note that the only allowed pole is at k 2 = 0, which has the meaning of approaching the boundary of 
the region Q. The partition function restricted to f l then becomes, 


Z(J) = / [dA] exp [—S F p] 


= / [dA] [dc] [dc] V(Cl) exp 


-Sym - j dx (c a d,D ab c b - ^(^A“) 2 


( 22 ) 


with 


nn) = e(i-a(0,A)), 


(23) 
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where 0(1 — cr(0, A)) is the Heaviside step function ensuring the no-pole condition. Considering the transversality 
of the gauge field in the Landau gauge and by making use of the integral representation of the Heaviside step 
function, one gets the following expression of the partition function restricted to the first Gribov region Q , 


Z{J) = A f [ / •DAe ^ 1 - <7(0 ’ A )) e -' Sp p , 

J 27Tld J 


with the ghost form factor given by 


, 1 1 Ng 2 

^A) = -- 


VdN 2 - 1 J (2tt) 


^ Ai(- q )A e a (q)± . 


(24) 


(25) 


At the end, the final expression for the gauge fixed Yang-Mills action accounting for the Gribov copies reads 


S' = S YM + S gi + S, 


(26) 


with Sp = /3 (cr(0 ,A) — 1). As mentioned before, our analytical work shall be restricted to first loop order and 
therefore we need only to consider quadratic terms of the quantized fields in the action given in equitation (26). 
Namely, we have 


with 


J quad 


•yab 




= <5 a6 + - - 1 + 


2Ng 2 /3 6^ 

(TV 2 - 1 )Vd q 2 


By performing the Gaussian integration of eq.(27) and making use of the relation 

(det V° 

one can finally get 


>“&) 2 _ e -|lndetP“J _ g —|Tr \nV* 


p+ioo+e 


J quad 


' —ioo+e 


c .f(0) 

2m 


where 


f(/3) = ^-hi/3- d -±(N 2 ~l)V / 111 (<? 2 + 


pNg 2 2 1 
TV 2 - 1 dVq 2 


(27) 


(28) 


(29) 


(30) 


(31) 


The factors (TV 2 — 1) and (d — 1) in front of the integral came from the trace over the SU(N) gauge group indices 
and from the trace over the Euclidean space-time indices, respectively 2 . Note that the factor (d— 1) is obtained only 
after the Landau gauge limit is taken. In the thermodynamic limit (when V —> oo) the saddle-point approximation 
becomes exact, so that the integral (30) can be easily computed, resulting in 


with /?* solving the equation 


«-V£„ 


df_ 

dp 


= e 


,/(/?*) 


= 0, 


(32) 


(33) 


which lead us to the gap equation 


d — 
d 



d d q 1 

(27T) d („4 , 2/3* Ng 2 \ 

[p "T" (N 2 -l)dV) 


1 . 


(34) 


Note that the Gribov parameter P, introduced to get rid of gauge ambiguities by restricting the path integral 
to the first Gribov region f2, is not a free parameter, it is given dynamically by solving its gap equation (34). 
Besides, it has dimension of [mass] 4 and is also proportional to the space-time volume V. Consequently, in the 
thermodynamic limit the logarithmic term of eq.(31) becomes null resulting in eq.(34). The function f{P*) can be 
interpreted, from eq.(32), as the free-energy of the system. 

x For more details concerning the derivation of the partition function restricted to the first Gribov region take a look at [ 0, 11, 12]. 

2 A careful computation of the functional trace of 1 n 'PY can be found in [12] 
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2.4 The gluon propagator 

In the present subsection we motivate that a possible sign of confinement could be read off from the poles of 
the gluon propagator, putting this quantity at the centre of any further discussion in the present work. At one-loop 
order only quadratic terms of the action (26) really matter, so that one can read off the two point function of the 
gauge field from the inverse of the operator (28), setting £ —► 0 at the very end of the computation. Namely, 


<A“( g )4(-g)) 


q 


2 


q 4 + 


2Nq 2 (3* 

(N 2 -l)dV 




Things become easier to analyze if a new Gribov’s parameter is defined, 


4 _ 2 (3*Ng 2 

7 (V 2 - 1 )dV ’ 


(35) 


(36) 


so that the gluon propagator can be decomposed, 


<A“(g)4(-g)) 


1 

2 


1 

q 2 + i y 2 




(37) 


Observe from (37) that the gluon propagator is suppressed in the infrared (IR) regime, while displaying two 
complex conjugate poles, m± = ±ky 2 . That feature does not allow us to attach the usual physical particle 
interpretation to the gluon propagator. From the analytic point of view the gluon propagator (37) has not a(n) 
(always) positive Kallen-Lehmann spectral representation, which is necessary to attach a probabilistic interpretation 
to the propagator 3 . These features lead us to interpret the gluon as being confined. 


3 The Yang-Mills + Higgs field theory 

Before properly starting to analyze the proposed model, let us state a few words on general features 
that are useful in this work. As mentioned before, the considered SU(2) and SU(2) x [7(1) Yang-Mills gauge 
theories are coupled to a scalar Higgs field. The Higgs field is considered in either the fundamental and the adjoint 
representation: in the SU( 2) case both the fundamental and the adjoint representation are analysed; while in the 
SU( 2) x C7(l) case only the the fundamental representation will be considered. 

Usually the unitary gauge arises as a good choice when the Higgs mechanism is being treated, since in 
this gauge physical excitations are evident. However, instead of fixing the unitary gauge, we are going to choose 
the more general gauge, whereby the unitary gauge is a special limit, £ —► oo, and the Landau gauge can be 
recovered when £ —> 0. In the case of interest, the Landau gauge is imposed at the end of each computation. 

The R^ gauge condition reads, 

f a = - !(^Vra(l)(T 0 ) m „V„ . (38) 

m,n 

Note that in eq. (38) iip m (x)(T a ) rnn v n is a possible form for the b a field, defined in (7). The <p field is defined as a 
small fluctuation around the vacuum configuration of the Higgs field $, 

$(x) = ip(x) + v , (39) 

with vacuum expectation value (ip) = 0. In order to apply the gauge condition (38) we should follow the standard 
process described in many text books, [26, 27, 28]. Rather than a Dirac delta function, as in eq. (4), one should 
put a Gaussian function, 


S(f a (A)) -> exp (-1 J d d x /“/“) ■ (40) 

In the limiting case of £ —> 0 the Gaussian term (40) oscillates very fast around /“ = 0 so that the Gaussian term 
(40) behaves like a delta function, ensuring the desired Landau gauge. On the other side, if £ —► oo then we have 

3 See [13, 14] and references therein for more details on the confinement interpretation of gluons, i-particles and the existence of 
local composite operators, related to these z-particles, displaying positive Kallen-Lehmann spectral representation. For lattice results 
pointing to the same confinement interpretation see [23, 24, 25]. 
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the unitary gauge. Needless to say, the limit £ —► 0, recovering the Landau gauge, should be applied at the very 
end of each computation 4 . 

In the present work we should deal with the Higgs field frozen at its vacuum configuration — i.e., $ = v. 
It is equivalent to replacing every Higgs field $ in the action with its vacuum expectation value v. Since the 
gauge fixing condition (38) only depends on the vector gauge field A M and on the fluctuation of the scalar field p 
multiplied by the gauge parameter £, the Gribov procedure remains valid for the limit case £ —> 0. 


4 The SU( 2)+Higgs field 

4.1 The Higgs field in the Fundamental representation: d-dimensional results 

In the present section nonperturbative effects of the 5t/(2)+Higgs model will be considered, by taking into 
account the existence of Gribov copies. The fundamental representation of the Higgs field, in d = 3 and d = 4, will 
be studied first. Subsequently, its adjoint representation, in d = 3 and d = 4, will be considered. 

Working in Euclidean spacetime, the starting action of the current model reads 

S = J d d x + X 2 ($t$ - v 2 f - {df ^l + ed.DfA , (41) 

where the covariant derivative is given by 

= d/F - > (42) 

and the vacuum expectation value of the scalar field is = v5 21 , with i = 1,2, so that all components of the 

gauge field acquire the same mass, m 2 =* a -^-- 

By following the procedure described in the subsection 2.2, the restriction to the first Gribov region S2 
relies on the computation of the ghost form factor, given by equation (25), and on the enforcement of the no-pole 
condition, (21)-(23). Since the presence of the scalar Higgs field does not influence the procedure of quantizing the 
gauge field, the non-local Gribov term, which is proportional to /?, is not affected and one ends up with the action 

S' = S + 0*a(O,A) -/?* , (43) 

with S given by eq.(41), the ghost form factor given by (25) and (3* stands for the Gribov parameter that solves 
the gap equation. 


4.1.1 The d-dimensional gluon propagator and the gap equation 

In order to compute the gluon propagator up to one-loop order in perturbation theory let us follow the 
steps described in the subsection 2.4. Therefore, we have 


(A;(q)Al(-q)) 


tab 




q 2 + 


4 g 2 /3* 

3 dV 




(44) 


whereby /?* solves the gap equation, which is obtained by the means of the subsection 2.3. After computing the 
Gaussian integral of the partition function and taking the trace over all indices, one ends up with the following 
partition function, 

Z quad = j ^ e^) = e~ v£v , (45) 

whereby one reads the free-energy, 


m 


/? — In/? 


3(d- 1)V 
2 


d d k 

(2n) d 


In 



4 g 2 /3 1 \ 

3 dV k 2 J ’ 


(46) 


4 lt is perhaps worthwhile pointing out here that the Landau gauge is also a special case of the ’t Hooft R £ gauges, which have 
proven their usefulness as being renormalizable and offering a way to get rid of the unwanted propagator mixing between (massive) 
gauge bosons and associated Goldstone modes, ~ A^d^(j). The latter terms indeed vanish upon using the gauge field transversality. 
The upshot of specifically using the Landau gauge is that it allows to take into account potential nonperturbative effects related to the 
gauge copy ambiguity. 












which is equivalent to (31). In the thermodynamic limit the integral of equation (45) may be solved through the 
saddle-point approximation (33) leading to the gap equation 5 


2(d-l) 2 f d d q 1 

5 J 


(47) 


Let us now carry out the analysis of the gluon propagator for d = 3 and d = 4, after solving the gap 
equation. 


4.1.2 The d = 3 case 

Now, let us proceed with the solution of the gap equation (47). Since we are working in d 
equation contain a finite integral, easy to be computed, leading to 

V _ 1 g 4 V 

3 dV P 4 \ 2 9t t 2 J ' 


3, the gap 

(48) 


As done in the previous section, the analysis of the gluon propagator could be simplified by making explicit use of 
its poles. Namely, 


with 


(A^AK-q)) 


6 ab 

9 O - 

— m_ 





i I'sV 2 



(49) 


(50) 


In this way, we may distinguish two regions in the (i/ 2 ,ff 2 ) plane: 

i) when g 2 < 9n 2 u 2 both masses (m 2 i _,m 2 _) are positive, as well as the residues. The gluon propagator, eq.(49), 
decomposes into two Yukawa modes. However, due to the relative minus sign in expression (49) only the 
heaviest mode with mass m 2 + represents a physical mode. We see thus that, for g 2 < 97r 2 z/ 2 , all components 
of the gauge field exhibit a physical massive mode with mass rr&_. This region is what can be called a Higgs 
phase. 

Let us also notice that, for the particular value g 2 = ^-v 2 , corresponding to a vanishing Gribov parameter 
/3 = 0 , the unphysical Yukawa mode in expression (49) disappears, as m 2 = 0. As a consequence, the gluon 
propagator reduces to that of a single physical mode with mass ^-u A . 

ii) when g 2 > 9n 2 is 2 , the masses (rri^,m‘L) become complex. In this region, the gluon propagator, eq.(49), 
becomes of the Gribov type, displaying complex conjugate poles. All components of the gauge field become 
thus unphysical. This region corresponds to the confining phase. 


4.1.3 The d = 4 case 

With quite the same process as for the d = 3 case, let us analyse the poles of the gauge field propagator 
by solving the gap equation for d = 4. To that end the following decomposition becomes useful 


! ^~ q2 + Y 13 = (l 2 + m\){q 2 +m 2 _) 


(51) 


with 


2 1 
m + = 2 


g 2 v 2 


5 4 A 


4 g 2 


m_ = - 


g 2 v 2 


g 4 A _ 

4 3 p 


(52) 


Making use of the MS renormalization scheme in d = 4 — e the gap equation (47) becomes 


In 


~P 2 


In 


P 2 


32?r 2 

w 


(53) 


5 We remind here that the derivative of the term ln/3 in expression (46) will be neglected, for the derivation of the gap equation, 
eq.(47), when taking the thermodynamic limit. 
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After a suitable manipulation we get a more concise expression for the gap equation 

2^/W ln(o) = - (l + In (l + v/w) + (l - V^?) In (l - v'w) , 


where we have introduced the dimensionless variables 


a = 


2 2 
9 v 




4/z 2 eU 3 » 


c = 


16 p* 

3 g 2 v A 


> 0 , 


(54) 


(55) 


with 0 < C, < 1 in order to have two real, positive, distinct roots {m 2 + ,m 2 _). For £ > 1, the roots (m+,rrvL) become 
complex conjugate, and the gap equation takes the form 

2 \J( — 1 ln(a) = —2 arctan (•*/£ — 1 ) — — 1 In ( . (56) 


Moreover, it is worth noticing that both expressions (54), (56) involve only one function, i.e. they can be written 
as 


2 In (a) = g{() , 


(57) 


where for g(£) we might take 


9(0 = (- (l + Vh) In (l + x/w) + (l - \/w) In (l ^ \/w)) , (58) 


which is a real function of the variable £ > 0. Expression (56) is easily obtained from (54) by rewriting it in the 
region £ > 1. In particular, it turns out that the function g(Q < —2 In 2 for all ( > 0, and strictly decreasing. As 
consequence, for each value of a < equation (57) has always a unique solution with ( > 0. Moreover, it is easy 
to check that g(l) = —2. Therefore, we can distinguish ultimately three regions, namely 

(a) when a > |, eq.(57) has no solution for Since the gap equation (47) has been obtained by applying the 
saddle point approximation in the thermodynamic limit, we are forced to set j3* = 0. This means that, when 
a > h, the dynamics of the system is such that the restriction to the Gribov region cannot be consistently 
implemented. As a consequence, the standard Higgs mechanism takes place, yielding three components of 
the gauge field with mass m 2 = Note that, for sufficiently weak coupling g 2 , a will unavoidably be 

larger than 


(b) when ^ < a < equation (57) has a solution for 0 < £ < 1. In this region, the roots (rn 2 + ,m 2 _) are real and 


the gluon propagator decomposes into the sum of two terms of the Yukawa type: 


(A%{q)A b v {- q )) = 


S ab 


m± — m_ \ q z + m 






(59) 


Moreover, due to the relative minus sign in eq.(59) only the component J-+ represents a physical mode. 


(c) for a < -, equation (57) has a solution for £ > 1. This scenario will always be realized if g 2 gets sufficiently 
large, i.e. at strong coupling. In this region the roots (mj_, rn _) become complex conjugate and the gauge 
boson propagator is of the Gribov type, displaying complex poles. As usual, this can be interpreted as the 
confining region. 


In summary, we clearly notice that at sufficiently weak coupling, the standard Higgs mechanism will definitely take 
place, as a > whereas for sufficiently strong coupling, we always end up in a confining phase because then a < 

Having obtained these results, it is instructive to go back where we originally started. For a fundamental 
Higgs, all gauge bosons acquire a mass that screens the propagator in the infrared. This effect, combined with a 
sufficiently small coupling constant, will lead to a severely suppressed ghost self energy, i.e. the average of (25) 
(to be understood after renormalization, of course). If the latter quantity will a priori not exceed the value of 1 
under certain conditions — i.e., satisfying the no-pole condition — the theory is already well inside the Gribov 
region and there is no need to implement the restriction. Actually, the failure of the Gribov restriction for a > ^ 
is exactly because it is simply not possible to enforce that cr(0) = 1. Perturbation theory in the Higgs sector is in 
se already consistent with the restriction within the 1st Gribov horizon. Let us verify this explicitly by taking the 
average of (25) with, as tree level input propagator, a transverse Yukawa gauge field with mass m 2 = s -^~. Using 
that there are 3 transverse directions 6 in 4 d, we have 

*(°) = 1 - ^2 ln ( 2fl ) ' ( 60 ) 


6 We have been a bit sloppy in this paper with the use of dimensional regularization. In principle, there are 3 — e transverse 
polarizations in d = 4 — e dimensions. Positive powers in e can (and will) combine with the divergences in e -1 to change the finite 
terms. However, as already pointed out before, a careful renormalization analysis of the Gribov restriction is possible, see e.g. [29, 12] 
and this will also reveal that the “1” in the Gribov gap equation will receive finite renormalizations, compatible with the finite 
renormalization in e.g. cr(0), basically absorbable in the definition of a. 
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For a > the logarithm is positive and it is then evident that cr(0) will not cross 1, indicating that the theory 
already is well within the first Gribov horizon. 

Another interesting remark is at place concerning the transition in terms of a varying value of a. If a crosses 
I, the imaginary part of the complex conjugate roots becomes smoothly zero, leaving us with 2 coinciding real 
roots, which then split when a grows. At a = one of the roots and its accompanying residue vanishes, to leave 
us with a single massive gauge boson. We thus observe that all these transitions are continuous, something which 
is in qualitative correspondence with the theoretical lattice predictions of the classic work [ ] for a fundamental 
Higgs field that is “frozen” (A —> oo). Concerning the somewhat strange intermediate phase, i.e. the one with a 
Yukawa propagator with a negative residue, eq.(59), we can investigate in future work in more detail the asymptotic 
spectrum based on the BRST tools developed in [30] when the local action formulation of the Gribov restriction 
is implemented. Recent works on the lattice confirm the existence of a cross-over region, where there is no line 
separating the “phases”, as e.g. [31, 32] where the authors work in the non-aligned minimal Landau gauge and 
observe the transition between a QCD-like phase and a Higgs-like phase, in a region away from the cross-over 
region. 


4.1.4 The vacuum energy, phase transition and how trustworthy are the results? 

Let us look at the vacuum energy £ v of the system, which can easily be read off from expression (45), 


namely 


£ v — —(3* + 


drk 

(2^ 


In k 2 + 


3 fc 2 


+ 


(61) 


where j3* is given by the gap equation (47). Making use of the MS renormalization scheme, the vacuum energy 
may be written as: 


for a < we have 
8 

f — 


l 


9 ff 4 i/ 4 


+ 


32t r 2 
1 1 
4 32tt 2 


1 - 


32t 


1 c 


i i 


2 32t r 2 4 32 tt 2 


(4-2£)(ln(a)-- 


(62) 


((l + V / l 3 C) 2 ln(l + x /w) + (l- v ^c) 2 ln(l- v ^c)^) , 


where £ is obtained through eqs.(57),(58). 
for a > 


9g 4 v 4 


f — 

C"?; 


32t r 2 


1 - 


32tC 


32t r 2 


ln («) - o + 


32t r 2 


In 2 . 


(63) 


From these expressions we could check that the vacuum energy £ v (a) is a continuous function of the variable a, as 
well as its first and second derivative, and that the third derivative develops a jump at a = 4. We might be tempted 
to interpret this is indicating a third order phase transition at a = The latter value actually corresponds to a line 
in the (g 2 , v) plane according to the functional relation (55). However, we should be cautious to blindly interpret 
this value. It is important to take a closer look at the validity of our results in the light of the made assumptions. 
More precisely, we implemented the restriction to the horizon in a first order approximation, which can only be 
meaningful if the effective coupling constant is sufficiently small, while simultaneously emerging logarithms should 
be controlled as well. In the absence of propagating matter, the expansion parameter is provided by y = 
as in pure gauge theory. The size of the logarithmic terms in the vacuum energy (that ultimately defines the 

gap equations) are set by In -=L and m 2 _ In-=7. A good choice for the renormalization scale would thus be 
Jl 2 ~ |m^|: for (positive) real masses, a fortiori we have m 2 _ < m 2 + and the second log will not get excessively large 
either because m_ gets small and the pre-factor is thus small, or m?_ is of the order of m 2 , and the log itself small. 
For complex conjugate masses, the size of the log is set by the (equal) modulus of m 2 ± and thus both small by our 
choice of scale. 

Let us now consider the trustworthiness, if any, of the a = \ phase transition point. For a ~ 4, we already 
know that ( ~ 0, so a perfect choice is ~p 2 ~ m\ ~ s -%~- Doing so, the a-equation corresponds to 


1 

2 


~ e 


-i+ 


4 

3 y 


(64) 


so that y ~ 4. Evidently, this number is thus far too big to associate any meaning to the “phase transition” at 
a = Notice that there is no problem for the a small and a large region. If v 2 is sufficiently large and we set 
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/ 1 2 ~ we have a small j/, leading to a large a, i.e. the weak coupling limit without Gribov parameter and 

normal Higgs-like physics. The logs are also well-tempered. For a small u 2 , the choice JI 2 ~ \Jg 2 9* will lead to 

1 I 4 

a ~ (small number)e~ + (65) 


so that a small a can now be compatible with a small y , leading to a Gribov parameter dominating the Higgs 

2 2 9 

induced mass, the “small number” corresponds to 9 A — . Due to the choice of Jr, the logs are again under control 


in this case. 

Within the current approximation, we are thus forced to conclude that only for sufficiently small or large 
values of the parameter a we can probe the theory in a controllable fashion. Nevertheless, this is sufficient to 
ensure the existence of a Higgs-like phase at large Higgs condensate, and a confinement-like region for small Higgs 
condensate. The intermediate a-region is more difficult to interpret due to the occurrence of large logs and/or 
effective coupling. Notice that this also might make the emergence of this double Yukawa phase at a = - ss 0.37 
not well established at this point. 


4.2 The Higgs field in the ajoint representation: (/-dimensional results 

Let us now assume the adjoint representation for the Higgs field in the action (41). In that case, the 
vacuum configuration which minimizes the energy is achieved by a constant scalar field satisfying 

($“) = vd a3 , (66) 

leading to the standard Higgs mechanism, so that the off-diagonal components of the gluon field A“, a = 1,2, 
acquires a mass m 2 H = g 2 v 2 , i.e. 


(A^p)A^-p)) 


5 a P 

p 2 + m 2 H 



P/j.pA 
P 2 ) 


while the third component A 3 ^ should remain massless, 


(67) 


(Alml(-P)) = £ (A- - ■ (68) 

However, as was pointed out by Polyakov [1], the theory exhibits a different behaviour. The action (41) 
admits classical solitonic solutions, known as the ’t Hooft-Polyakov monopoles' which play a relevant role in the 
dynamics of the model. In fact, it turns out that these configurations give rise to a monopole condensation at weak 
coupling, leading to a confinement of the third component A 3 , rather than to a Higgs type behaviour, eq.(68), a 
feature also confirmed by lattice numerical simulations [5, 6]. 

Since our aim is that of analysing the nonperturbative dynamics of the Georgi-Glashow model by taking 
into account the Gribov copies, let’s follow the procedure described in the subsection 2.3. Due to the presence of 
the Higgs field in the adjoint representation, causing a breaking of the global gauge symmetry, the ghost two-point 
function has to be decomposed into two sectors, diagonal and off-diagonal: 


where 


g J k; 


A) 


Soff(k-A) - fc2 (l + cjoff{k-A))^ fc2 > 


Qdiagik] A) - k2 (l + a dlag (k;A))^ fc2 ( x A) ) 


(69) 


(70) 


(71) 


As we know, the quantities a 0 ff(k;A), <Jdiag{k\A ) turn out to be decreasing functions of the momentum 
k and making use of the gauge field transversality, we have 


a °ff( 0; A) 
e^diag (0> A) 


f d d q (Al(q)Al(-q) + \A°(q)A“(-q)) 
Vd J (2Tr) d q 2 

g 2 f d d q (A“(q)A“(-g)) 

VdJ (2n) d q 2 


7 These configurations are instantons in Euclidean space-time. 


(72) 
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Therefore, the no-pole condition for the ghost function Q ab (k,A) is implemented by imposing that [LO, 12, 11] 


<Toff(0\ A) < 1 , 

&diag(Q'i A) 1 . 


(73) 


After that we need two different parameters to implement the no-pole condition in the action, leading us 
to the following action accounting for the Gribov ambiguities, 


S ' = S + /3* (a off ( 0, A) - 1) + w* {a diag (0, A) - 1) . 
In the action (74) /3* and uj* are given dynamically through its own gap equation. 


(74) 


4.2.1 The (-/-dimensional gluon propagator and the gap equation 

The partition function accounting only for quadratic terms of the action (74) can be written as 


A, 


d/3e^dwe“ rj „ airj <3l ^ 


2ni/3 2 ttiuj 


— [dA Q ][dA 8 ] e 2 J <*« 


(75) 


with 


KS = <” s (v(9 2 + -V) + ()-i)^ + ^(^ + |)4i 


HIS I , 


s 


■[AV 


= duv q - 


2 2 Ujg 2 1 


Vd q 2 


+ 


- 1 


9m 


(76) 


The parameter £ stands for the usual gauge fixing parameter, to be put to zero at the end in order to recover the 
Landau gauge. Evaluating the inverse of the expressions (76) and taking the limit £ —► 0, the gluon propagators 
become 


(Al(q)Al(-q)) = 


4 

9 ^ Vd 


S,,u - 


q^ 


<A“(g)A£(-g)> = 5^ 


9 2 (q 2 + g 2 " 2 ) + 2 ih (P + l) 


A q» 

OfAv o 

<r 


(77) 

(78) 


The off-diagonal sector of the gluon propagator can be put in a more convenient form by making explicit its poles. 
Namely, 


with 


(A^q)A^-q)) = 


gW + y/fiA-lr 


5 al3 


m+ — m_ \q + q 2 + m_ 


- 


9m q» 


2 9 

m_ = 


2 i/ 2 - \/g 4 v 4 - 4r 


r = 


2 g 2 




UJ 

+ 2 


(79) 


(80) 


2 ’ 2 Vd 

Since the Gribov parameters (/?, uj) are fixed dynamically through the gap equation, now we should integreat 
out the gauge field from equation (75) and make use of the saddle-point approximation, in the thermodynamic 
limit, which will gives us two gap equations, enabling us to express /3 and ui in terms of the parameters of the 
starting model, i.e. the gauge coupling constant g and the vev of the Higgs field u. Performing the Gaussian 
integral of the partition function (75) we get 8 


Zquad 


_ [ *) 


27 xi 2iri 


(81) 


with 


f(ui,/3) = (3 + uj — ln/3 — lnw — 


(d- 1)V f d d q 


( 2n)° 


In q z + 


2 2uig 2 1 


2(d-l)V r d d q 

~2 J j2nY 


In (^q 2 + g 2 v 2 ) + ^(2/3 + 0 ,)^ 


Vd q 2 

1 


(82) 


8 In order to obtain eq. (82) we had to take the functional trace of \nV and In Q following the process described in subsection 2.3. 
For more details take a look at [12]. 


13 

























Since in the thermodynamic limit, as mentioned in the section 2.3, the integral (81) can be solved through the 
saddle point approximation, 


df 

dp* 



(83) 


one gets the following two gap equations 


4(d — l)g 2 r d d q 1 

2d J (2n) d g 4 + 2 wfZ 


4(d-l)<? 2 

f d d q j 

< 

2d J 

(27r) d ! 

\q 2 (q 2 +g 2 v' 2 ) + g 2 1 

(T + S) 



(84) 

(85) 


Therefore, p* and uj* can be expressed in terms of the parameters v,g. To solve the gap equations the denominator 
of eq.(85) can be decomposed into its poles, which is similar to (78)-(80). 


4.3 The d — 3 case 

In the three-dimensional case both gap equations cause not many difficulties to be solved, as there are no 
divergences to be treated. Namely, the first gap equation, eq.(84), leads to the following result, 


w*(ff) = 

while the second one, given by eq.(85), leads to 


T = 


,2g 2 * g 2 


9 


1 2 2 ff 4 

-q v -- 

2 y 32 tt 2 


( 86 ) 


(87) 


Now we can look at the gluon propagators, (77) and (78), and analyse the different regions in the (g,v) 
plane. Let us start by the diagonal component A 3 . Namely, we have 


( A l(q) A l(-q)) 


q 4 + 






( 88 ) 


One observes that expression ( 88 ) turns out to be independent from the vev v of the Higgs field, while displaying 
two complex conjugate poles. This gauge component is thus of the Gribov type. In other words, the mode A ^ is 
always confined, for all values of the parameters g, u. Concerning now the off-diagonal gluon propagator (78), after 
decomposing it into two Yukawa modes (79), we could find the following regions in the ( g , v) plane: 


i) when g 2 < 327 r 2 ^ 2 , corresponding to r < 3 ” , both masses to^_, m 2 _ are real, positive and different from each 
other. Moreover, due to the presence of the relative minus sign in expression (79), only the heaviest mode 
with mass m 2 + represents a physical excitation — i.e ., despite the existence of two real positive poles, m 
and wP _, only the contribution related to the m 2 , pole has physical meaning. 

It is also worth observing that, for the particular value g = 16tt 2 i/ 2 , corresponding to r = 0, the unphysical 
mode in the decomposition (79) disappears. Thus, for that particular value of the gauge coupling, the 
off-diagonal propagator reduces to a single physical Yukawa mode with mass 167 t 2 zA 

ii) when g 2 > 32tP 1 v 2 , corresponding to r > all masses become complex and the off-diagonal propagator 

becomes of the Gribov type with two complex conjugate poles. This region, called Gribov region since all 
modes are of Gribov type, corresponds to a phase in which all gauge modes are said to be confined. 


In summary, when the Higgs field is in the adjoint representation we could find two distinct regions. For g 2 < 327 t 2 v 2 
the H 3 mode is confined while the off-diagonal propagator displays a physical Yukawa mode with mass m^_. This 
phase is referred to as the U( 1) symmetric phase [5, 6 ]. When g 2 > 32n 2 v 2 all propagators are of the Gribov 
type, displaying complex conjugate poles leading to a confinement interpretation. According to [5, 6 ] this regime 
is referred to as the SU{2) confined phase. 

Since our results were obtained in a semi-classical approximation (i.e., lowest order in the loop expansion), 

let us comment on the validity of such approximation. In general, the perturbation theory is reliable when the 

2 

effective coupling constant is sufficiently small. The effective coupling depends, in 3d, on the factor ^ 4 ^ 3/2 • However, 
since g 2 has mass dimension 1 the effective coupling is not complete yet. In the presence of a mass scale M, the 
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perturbative series — for e.g. the gap equation — will organize itself automatically in a series in G 2 /M. Let 
us analyse, for example, the case where g 2 < 327 t 2 v 2 , the called the “Higgs phase”. In this case the effective 
coupling will be sufficiently small when is small compared 9 to 1. Such condition is not at odds with the 

retrieved condition g 2 < 327 t 2 v 2 . Next, assuming the coupling g 2 to get large compared to v 2 , thereby entering 
the confinement phase with cc masses, g 2 dominates everything, leading to a Gribov mass scale r oc g 8 , and an 
appropriate power of the latter will secure a small effective expansion parameter consistent with the condition 
g 2 > 32n 2 is 2 . We thus find that at sufficiently small and large values of ^ our approximation and results are 
trustworthy. 


4.4 The d = 4 case 

Let us start by considering the second gap equation, eq.(85). Performing the decomposition described in 
eq.(80) the referred gap equation becomes of the same form as the one obtained in the fundamental d = 4 case, 
eq.(53). The difference between the fundamental and adjoint d = 4 cases appears in the definition of the mass 
parameter rn± (see eq.(52) and eq.(80), respectively). Namely, 


1 + 


In 


mi — m_ 


n 2 


Introducing now the dimensionless variables 10 


fo = 


g 2 v 2 


1 


mi — m_ 


g 2 v 2 


In 


and 


g 2 


32t r 2 
3 g 2 


£ = 


4 T 

7^ i0 - 


(1 32tt 2 X (. 272-tt 2 A /\2 ’ 

2 fi 2 e'' 3 s 2 / 2 ev 21 s 2 / MS 

with 0 < £ < 1 . Proceeding as in the fundamental d = 4 case, eq. (89) can be recast in the following form 


(89) 


(90) 


2 V / W ln(fo) = - (l + x/w) In (l + x/w) + (l - x/w) In (l - x/w) ■ (91) 


or compactly, 


2Info = 5(0 • 


(92) 


Also here, in the adjoint d = 4 case, eq.(92) remains valid also for complex conjugate roots, viz. £ > 1. We are 
then led to the following cases. 


4.4.1 When fo < 7 


Using the properties of g(£), it turns out that eq.(91) admits a unique solution for £, which can be explicitly 
constructed with a numerical approach. More precisely, when the mass scale g 2 v 2 is sufficiently smaller than AjU-, 

i.e. 


2 2 

g » 


< 2 eV 


(i- 3 ©?) A— 


(93) 


we have what can be called the 17(1) confined phase. In fact, in this regime the gap equation (84) leads to a non¬ 
null ui*, so that the diagonal component of the gauge field is said to be of the Gribov type, i.e. with confinement 
interpretation. 

On the other side, the second gap equation (91) splits this region in the two subregions: 

(i) when 7 < fo < \ equation (91) has a single solution with 0 < £ < 1. In this region, the roots (m^, m 2 _) are 
thus real and the off-diagonal propagator decomposes into the sum of two Yukawa propagators. 

However, due to the relative minus sign in eq.(79), only the component with m 2 , pole can be associated to a 
physical mode, analogously as in the fundamental case. Due to the confinement of the third component A 3 , 
this phase is recognized as the 17(1) confining phase. It is worth observing that it is also present in the 3 d 
case, with terminology coined in [5], see also [18]. 

(ii) for fo < 7, equation (91) has a solution for £ > 1. In this region the roots (m^_, m 2 _) become complex conjugate 
and the off-diagonal gluon propagator is of the Gribov type, displaying complex poles. In this region all gauge 
fields display a propagator of the Gribov type. This is recognized as the SU{2) confining region. 


Similarly, the above regions are continuously connected when fo varies. In particular, for fo 7, W e obtain £ = 0 as 
solution. 

9 The Higgs mass v 1 is then the only mass scale entering the game. 

10 We introduced the renormalization group invariant scale A Ms . 
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4.4.2 The case b > \ 


Let us consider now the case in which b > Here, there is no solution of the equation (91) for the 
parameter £, as it follows by observing that the left hand side of eq.(91) is always positive, while the right hand 
side is always negative. This has a deep physical consequence. It means that for a Higgs mass m 1 2 Higgs = g 2 v 2 
sufficiently larger than A^, i.e. 


2 2 
9 v 


>2 


i-- 


A| 


(94) 


the gap equation (85) is inconsistent. It is then important to realize that this is actually the gap equation obtained 
by acting with on the vacuum energy E v = — /(w, /3). So, we are forced to set j3* = 0, and confront the remaining 
w-equation, viz. eq.(84), which can be transformed into 


4 ln( 6 ) = 


(l + In (l + + (l - AA) In (l - 


- n/T 7 ! ln£ - V / Wln 2 l = ft(£) 


(95) 


after a little algebra, where £ = A?- The behaviour of /i(£) for £ > 0 is more complicated than that of c/(£). 
Because of the — ln£ contribution, /i(£) becomes more and more positive when £ approaches zero. In fact, /i(£) 
strictly decreases from +oo to — oo for £ ranging from 0 to +oo. 

It is interesting to consider first the limiting case b -4- yielding £ ss 1.0612. So, there is a discontinuous 
jump in £ (i.e. the Gribov parameter for fixed v) when the parameter b crosses the boundary value \. 

We were able to separate the b > \ region as follows: 


(a) For \ < b < 


0.51, we have a unique solution £ > 1, i.e. we are in the confining region again, with all 


2 y/vTe 

gauge bosons displaying a Gribov type of propagator with complex conjugate poles. 


(b) For 


V V^e 


< b < oo, we have a unique solution £ < 1 , indicating again a combination of two Yukawa modes 


for the off-diagonal gauge bosons. The “photon” is still of the Gribov type, thus confined. 


Completely analogous as in the fundamental case, it can be checked by addressing the averages of the expressions 
(72) that for b > | and lo obeying the gap equation with /? = 0, we are already within the Gribov horizon, making 
the introduction of the second Gribov parameter /3 obsolete. 

It is obvious that the transitions in the adjoint case are far more intricate than in the earlier studied 
fundamental case. First of all, we notice that the “photon” (diagonal gauge boson) is confined according to its 
Gribov propagator. There is never a Coulomb phase for b < oo. The latter finding can be understood again from 
the viewpoint of the ghost self-energy. If the diagonal gluon would remain Coulomb (massless), the off-diagonal 
ghost self-energy, cfr. eq.(72), will contain an untamed infrared contribution from this massless photon 11 , leading 
to an off-diagonal ghost self-energy that will cross the value 1 at a momentum k 2 > 0 , indicative of trespassing the 
first Gribov horizon. This crossing will not be prevented at any finite value of the Higgs condensate v, thus we are 
forced to impose at any time a nonvanishing Gribov parameter to. Treating the gauge copy problem for the adjoint 
Higgs sector will screen (rather confine) the a priori massless “photon”. 

An interesting limiting case is that of infinite Higgs condensate, also considered in the lattice study of 
13]. Assuming v — > oo, we have b — > oo according to its definition (90). Expanding the gap equation (95) around 
£ = 0 + , we find the limiting equation b 4 = |, or equivalently w* a h^ KS /g 4 v 4 . Said otherwise, we find that also the 
second Gribov parameter vanishes in the limit of infinite Higgs condensate. As a consequence, the photon becomes 
truly massless in this limit. This result provides —in our opinion— a kind of continuum version of the existence of 
the Coulomb phase in the same limit as in the lattice version of the model probed in [33]. It is instructive to link 
this back to the off-diagonal no pole function, see Eq. (72), as we have argued in the proceeding paragraph that 
the massless photon leads to <J 0 ff( 0) > 1 upon taking averages. However, there is an intricate combination of the 
limits v —» oo, w* —> 0 preventing such a problem here. Indeed, we find in these limits, again using dimensional 
regularization in the MS scheme, that 


cr o//( 0 ) 



128t r 2 

3 ff 2 

128tt 2 


A 1 f d 4 q 1 

(2A q 4 + + J (Ag 2 (g 2 + gV) + ^ 



WJT 

2 ~jJ 4 


In 


A 2 

n 2 




♦ 6 4 

lj g^v 

>778 


2 /T 




3ff 2 

647T 2 


ln8<? 2 + 


(96) 


1 The “photon” indeed keeps it coupling to the charged (= off-diagonal) ghosts, as can be read off directly from the Faddeev-Popov 

term c a di_ l D“ b c b . 
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The latter quantity is always smaller than 1 for g 2 positive, meaning that we did not cross the Gribov horizon. 
This observation confirm in an explicit way the intuitive reasoning also found in section 3.4 of [34], at least in the 
limit v —> oo. The subtle point in the above analysis is that it is not allowed to naively throw away the 2nd integral 
in the first line of (96) for v —> oo. There is a logarithmic \nv (v —» oo) divergence that conspires with the lnw* 
(w* —> 0) divergence of the 1st integral to yield the final reported result. This displays that, as usual, certain care 
is needed when taking infinite mass limits in Feynman integrals. 

4.4.3 The vacuum energy in the adjoint case 

As done in the case of the fundamental representation, let us work out the expression of the vacuum energy 
£ v , for which we have the one loop integral representation given by eq.(82) multiplied by —1. Making use of the 
MS renormalization scheme in d = 4e the vacuum energy becomes 


3£' 


g 4 v 4 
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"<?2 ~ 128tt 2 
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(ln(2& 2 0 - 1) 
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128tt 2 


(^(l^\/r^) 2 ln(l-+ (l + yi^) 2 ln(l + v ^^)^ , 


where b was introduced via its definition (90), while 


e' = 


At' 
g 4 v 4 1 


£ = 


At 

g 4 v 4 


and 


g 2 u* 



(97) 


(98) 


Since we found scenarios completely different for b < 1/2 and b > 1/2 with the scalar Higgs field in its 
adjoint representation, it becomes of great importance analysing the plot of the vacuum energy as a function of b. 
From Fig.3 one can easily find out a clear jump for b = 1/2, which can be seen as a reflection of the discontinuity 
of the parameter £. 



Figure 3: Plot of the vacuum energy in the adjoint representation as a function of the parameter b. The discontinuity at 
b = 1 is evident. 

Investigating the functional (97) in terms of £ and it is numerically (graphically) rapidly established 
there is always a solution to the gap equations ^ ^ = 0 for b < 4, but the solution £* is pushed towards the 

boundary £ = 0 if b approaches 4, to subsequently disappear for b > \ 12 . In that case, we are forced to return on 
our steps as in the fundamental case and conclude that /3 = 0, leaving us with a single variable £ = £' and a new 
vacuum functional to extremize. There is, a priori, no reason for these 2 intrinsically different vacuum functionals 
be smoothly joined at b = \. This situation is clearly different from what happens when a potential has e.g. 2 
different local minima with different energy, where at a first order transition the two minima both become global 
minima, thereafter changing their role of local vs. global. Evidently, the vacuum energy does not jump since it is 
by definition equal at the transition. 

Nevertheless, a completely analogous analysis as for the fundamental case will learn that b = \ is beyond 
the range of validity of our approximation 13 . The small and large b results can again be shown to be valid, so at 

2T he gap solutions correspond to a local maximum, as identified by analysing the Hessian matrix of 2nd derivatives. 

13 A little more care is needed as the appearance of two Gribov scales complicate the log structure. However, for small b the Gribov 
masses will dominate over the Higgs condensate and we can take JL of the order of the Gribov masses to control the logs and get a 
small coupling. For large b, we have /3* = 0 and a small ui*\ the first log will be kept small by its pre-factor and the other logs can be 
managed by taking JL of the order of the Higgs condensate. 
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large b (~ large Higgs condensate) we have a mixture of off-diagonal Yukawa and confined diagonal modes and at 
small b (~ small Higgs condensate) we are in a confined phase. In any case we have that the diagonal gauge boson 
is not Coulomb-like, its infrared behaviour is suppressed as it feels the presence of the Gribov horizon. 


5 The Electroweak theory: 577(2) x ?7(1)+Higgs field 


From now on in this work only the fundamental case of the Higgs field will be treated, for reasons relying 
on the physical relevance of the fundamental representation of this field. As a first step, we are going to present, 
as in the previous sections, general results for d-dimension. Afterwards, the 3 and 4-dimensional cases will be 
considered in the subsections 5.1 and 5.2. The starting action of the SU( 2) x [/(1)+Higgs field reads 


S = 


J? a JT' a 
4 M*' 




lAvBfii/ 
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=a 2 - ( d M 2 


2 £ 
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(99) 


where the covariant derivative is defined by 


= d^ 1 - - ig¥¥-Al& . ( 100 ) 

and the vacuum expectation value (vev) of the Higgs field is (<&*) = vS 2 ’. The indices i, j = 1,2 refer to the 
fundamental representation of SU( 2) and r“,a = 1,2,3, are the Pauli matrices. The coupling constants g and g' 
refer to the groups SU{2) and Z7(l), respectively. The field strengths and B^ are given by 

F% = d,A“ - d v Al + ge abc A\Al , B, lv = d\B V - d v B, t . (101) 

In order to obtain the boson propagators only quadratic terms of the starting action are required and, due 
to the new covariant derivative, this quadratic action is not diagonal any more. To diagonalize this action one 
could introduce a set of new fields, related to the standard ones by 

w+ = ^(K + iA l) > w 7 = ^ « - * A l) > ( 102a ) 

= y g 2 + g/2 (~ 9A ^ + 9 ' B v) and A » = y g2 + g/2 ( 9 ' A M + ■ ( 102 b) 

The inverse relation can be easily obtained. With this new set of fields the quadratic part of the action reads, 

S q uad = J d 3 x (\(d,w+ - d v w+)(d,w- - d v w~) + 

+ j d 3 x Q(0 M Z„ - d v Zrf + (g2 + p u2 Z ll Z„ + ^(d,A„ - d„A,) 2 ^J , (103) 


from which we can read off the masses of the fields W + , W , and Z: 


2 9 2 v 2 

m w = 


m 2 z = 


(g +9 y 


(104) 


The restriction to the Gribov region H still is needed and the procedure here becomes quite similar to what 
was carried out in the subsection 4.2. Due to the breaking of the global gauge invariance, caused by the Higgs 
field (through the covariant derivatives), the ghost sector can be split up in two different sectors, diagonal and 
off-diagonal. Namely, the ghost propagator reads, 


Q ab {k-,A) = 

By exploiting the ghost form factor we have 

Q 0 ff{k]A) 


5 a Pg off (k;A ) 0 

o Gdiag(k;A) 


k 2 \1 - a 0 ff(k; A) 


(105) 


(106) 
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and 


g diag (k-A) ~ k2 ( 1 _ a J agik . A) ) > 

where 

v°tt(0; A ) = ^ j^s ^ + A lM A V-pi) 

and 

a dtag(0', A) = ^a;(p)a;(- p ) . 

In order to obtain expressions (108a) and (108b), where V denotes the (infmte) space-time volume, the transversality 
of the gluon field and the property that a(k; A) a ff and cr(fc; A)di ga are decreasing functions of k were used 14 . From 
equations (106) and (108a) one can easily read off the two no-pole conditions. Namely, 

tr o //(0; A) < 1 , (109a) 

and 

CTdiag(0; A) < 1 . (109b) 

At the end, the partition function restricted to the first Gribov region Q reads, 

Z = f-^--^-[dA\[dB] . ( 110 ) 

J 27 riuj 27 rip 

5.0.4 The d-dimensional gluon propagator and the gap equation 

The perturbative computation at the semi-classical level requires only quadratic terms of the full action, 
defined in eq.(110) (with S given by eq.(99)), yielding a Gaussian integral over the fields. Inserting external fields 
to obtain the boson propagators, one gets, after taking the limit £ —»■ 0, the following propagators, 
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Moving to the fields WA,W^ ^Z^A^, one obtains 
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4 For more details concerning the ghost computation see [15, 16, 17, 12] 
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(111b) 

(111c) 

(Hid) 
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As expected, all propagators get deeply modified in the IR by the presence of the Gribov parameters /3 and oj. 
Notice, in particular, that due to the parameter oj a mixing between the fields A^ and Z M arises, eq.(112d). As such, 
the original photon and the boson Z loose their distinct particle interpretation. Moreover, it is straightforward to 
check that in the limit /3 —>■ 0 and w —>■ 0, the standards propagators are recovered. 

Let us now proceed by deriving the gap equations which will enable us to (dynamically) fix the Gribov 
parameters, (3 and oj, as function of g , g' and v 2 . Thus, performing the path integral of eq.(llO), in the semi-classical 
level, we get 


/(<*>,/?) 



2 (d - 1) f d d p 

2 J (2t rp 
1) 
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dV p 2 


log [p + y g 
J logA+(p,w)A_(p,w) 


(113) 


In eq.(113), /(w,/3) is defined according to eq.(81) and 


A± 


( p 4 + ^P 2 (a 2 + g' 2 ) + #) ± \J[^(g 2 +g' 2 )P 2 + 


dV 


^ v 2 g 2 g' 2 p 2 


(114) 


Making use of the thermodynamic limit, where the saddle point approximation takes place, we have the two gap 
equations given by 15 


and 


4(d-l) 2 f d d p 1 

2d 9 J (2n) d p 4 + a^l p 2 + W 

2(d-l) 2 f d d p _ p 2 + 4 ff' 2 _ 

d 9 J (2n) d p6 + ^ {g 2 + g/ 2 )pi + 2^2 + 


(115) 


(116) 


Given the difficulties in solving the gap equations (115) and (116), we propose an alternative approach 
to probe the gluon propagators in the parameter space v , g and g'. Instead of explicitly solve the gap equations, 
let us search for the necessity to implement the Gribov restriction. For that we mean to compute {u 0 ff( 0)) 
and (<Jdiag( 0)) with the gauge field propagators unchanged by the Gribov terms, be., before applying the Gribov 
restriction. Therefore, if (cr o //(0; A)) < 1 and (crdiagiO', A)) < 1 already in this case (without Gribov restrictions), 
then we would say that there is no need to restrict the domain of integration to f2. In that case we have, immediately, 
/3* = oj* = 0 and the standard Higgs procedure takes place. Namely, the expression of each ghost form factor is 


and 


off (0)) 


(d — l)g 2 r d d p 1 ( 1 

d J (2ir) d p 2 + yl g 2 
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p 2 + k r(g 2 + g 12 ) 



(o'diagi 0)) 


2(d — l)g 2 r d d p 1 / 1 \ 

d J {2n) d p 2 ^p2 + ylg 2 J 


(117) 


(118) 


5.1 The d = 3 case 

In the three-dimensional case things become easier since there is no divergences to treat. Therefore, 
computing the ghost form factors (117) and (118) we led to the following conditions 

(1 + cos(dw))~ < 3\/27t (119a) 

2^ < 3V2tt , (119b) 

where 9{W) stands for the Weinberg angle, 

C0S(M = -r ==,f ■ ( 12 °) 

vV + g ' 2 

These two conditions make phase space fall apart in three regions, as depicted in Figure 4. 

1 "'For more details see [ 5, 16, 17]. 
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Figure 4: There appear to be four regions in phase space. The region I is defined by condition (119b) and is characterized 
by ordinary Yang-Mills-Higgs behaviour (massive W and Z bosons, massless photon). The region II is defined 
by (119a) while excluding all points of region I — this region only has electrically neutral excitations, as the W 
bosons are confined (see Section 5.1.1); the massive Z and the massless photon are unmodified from ordinary 
Yang-Mills-Higgs behaviour. Region III has confined W bosons, while both photon and Z particles are massive 
due to influence from the Gribov horizon; furthermore there is a negative-norm state. In region IV all SU(2) 
bosons are confined and only a massive photon is left. Mark that the tip of region III is hard to deal with 
numerically — the discontinuity shown in the diagram is probably an artefact due to this difficulty. Details are 
collected in Section 5.1.2. 


• If g/v < 37 r/\/2, neither Gribov parameter is necessary to make the integration cut off at the Gribov horizon. 
In this regime the theory is unmodified from the usual perturbative electroweak theory. 

• In the intermediate case 3tt/\/ 2 < g/v < 3v / 27t/( 1 + cos Qw) only one of the two Gribov parameters, f3, is 
necessary. The off-diagonal {W) gauge bosons will see their propagators modified due to the presence of a 
non-zero /3, while the Z boson and the photon A remain untouched. 

• In the third phase, when g/v > 3\/2tt/{1 + cos0vy), both Gribov parameters are needed, and all propagators 
are influenced by them. The off-diagonal gauge bosons are confined. The behaviour of the diagonal gauge 
bosons depends on the values of the couplings, and the third phase falls apart into two parts, as detailed in 
section 5.2.2. 

Note that here in the 3-dimensional SU(2) x C/(1)+Higgs case, as well as in the 3 d >S'[/(2)-|-Higgs treated in section 
4.3, an effective coupling constant becomes of utmost importance when discussing the trustworthiness of the our 
semi-classical results. 


5.1.1 The off-diagonal {W) gauge bosons 

Let us first look at the behaviour of the off-diagonal bosons under the influence of the Gribov horizon. 
The propagator (112a) only contains the /3 Gribov parameter, meaning that to need not be considered here. 

In the regime g/v < 3tt/\/2 (region I in Figure 4) the parameter (3 is not necessary, due to the ghost form 
factor (<r diag (0)) always being smaller than one. In this case, the off-diagonal boson propagator is simply of massive 
type, with mass parameter \g 2 . 

In the case that g/v > 3tt/\/2 (regions II, III, and IV in Figure 4), the relevant ghost form factor is not 
automatically smaller than one any more, and the Gribov parameter (3 becomes necessary. The value of (3* is 
determined from the gap equations (115). After rewriting the integrand in partial fractions, the integral in the 
equation becomes of standard type, and we readily find the solution 


32 



( 121 ) 


Mark that, in order to find this result, we had to take the square of both sides of the equation twice. One can 
easily verify that, in the region g/v > 37 t/\/2 which concerns us, no spurious solutions were introduced when doing 


Replacing this value of f3* in the off-diagonal propagator (112a) one can immediately check that it clearly 
displays two complex conjugate poles. As such, the off-diagonal propagator cannot describe a physical excitation 
of the physical spectrum, being adequate for a confining phase. This means that the off-diagonal components of 
the gauge field are confined in the region g/v > 3t:/\/2. 
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5.1.2 The diagonal SU(2) boson and the photon field 

The other two gauge bosons — the A ^ and the B M — have their propagators given by (111b), (111c), and 
(llld) or equivalently — the Z M and the A^ - by (112b), (112c) and (112d). Here, oj is the only one of the two 
Gribov parameters present. 

In the regime g/v < 3v / 27r/(l + cos<V) (regions I and II) this uj is not necessary to restrict the region of 
integration to within the first Gribov horizon. Due to this, the propagators are unmodified in comparison to the 
perturbative case. 

In the region g/v > 3v / 2 , 7r/(l + cos 9w) (regions III and IV) the Gribov parameter w does become necessary, 
and it has to be computed by solving its gap equation, eq. (116). Due to its complexity it seems impossible 
to do so analytically. Therefore we turn to numerical methods. Using Mathematica the gap equation can be 
straightforwardly solved for a list of values of the couplings. Then we determine the values where the propagators 
have poles. 

The denominators of the propagators are a polynomial which is of third order in p 2 . There are two cases: 
there is a small region in parameter space where the polynomial has three real roots, and for all other values of 
the couplings there are one real and two complex conjugate roots. In Figure 4 these zones are labelled III and IV 
respectively. Let us analyze each region separately. 

5.1.3 Three real roots (region III) 

Region III is defined by the polynomial in the denominators of (111b), (111c), and (Hid) having three 
real roots. This region is sketched in Figure 4. (Mark that the tip of the region is distorted due to the difficulty in 
accessing this part numerically.) 

The residues of related to these poles were computed numerically. Only the two of the three roots have 
positive residue and can correspond to physical states. Those are the one with highest and the one with lowest 
mass squared. The third of the roots, the one of intermediate value, has negative residue and thus belongs to some 
negative-norm state, which cannot be physical. 

All three states have non-zero mass for non-zero values of the electromagnetic coupling g ', with the lightest 
of the states becoming massless in the limit g' —> 0. In this limit we recover the behaviour found in this regime 
in the pure SU{2) case [18] (the Z-boson field having one physical and one negative-norm pole in the propagator) 
with a massless fermion decoupled from the non-Abelian sector. 

5.1.4 One real root (region IV) 

In the remaining part of parameter space, there is only one state with real mass-squared. The two other 
roots of the polynomial in the denominators of (111b), (111c), and (llld) have non-zero imaginary part and are 
complex conjugate to each other. In order to determine whether the pole coming from the real root corresponds to 
a physical particle excitation, we computed its residue, which can be read off in the partial fraction decomposition 
(the result is plotted in Figure 5). It turns out the residue is always positive, meaning that this excitation has 
positive norm and can thus be interpreted as a physical, massive particle. The poles coming from the complex 
roots cannot, of course, correspond to such physical excitations. 

In the limit g' —> 0 we once more recover the corresponding results already found in the pure SU{2) case 
[ 8] (two complex conjugate poles in the propagator of the non-Abelian boson field) plus a massless photon not 
influenced by the non-Abelian sector. 

We shall emphasise here the complexity of the found “phase spectrum” in the 3d case. For the most part 
of the (g r /v, g/v) plane we found the diagonal component of the bosonic field displaying a mix of physical and non¬ 
physical contributions, regarding the regions III and IV. The off-diagonal component was found to have physical 
meaning only in the region I. The transition between those regions was found to be continuous with respect to the 
effective perturbative parameter ~ g/v. 

5.2 The d — 4 case 

In the 4-dimensional case the diagonal and off-diagonal ghost form factors read, using the standard MS 
renormalization procedure, 


/ , ,, 3 g 2 , 2a 
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Figure 5: The residue of the pole of the photon propagator. It turns out to be positive for all values of the couplings within 
the region IV. 


and Qw stands for the Weinberg angle. With expression (122) we are able to identify three possible regions, depicted 
in Figure 6: 

• Region I, where (crdiag(O)) < 1 and (cr o ff(0)) < 1, meaning 2a > 1. In this case the Gribov parameters are 
both zero so that we have the massive W^ 1 and Z , and a massless photon. That region can be identified with 
the “Higgs phase”. 

• Region II, where (<Tdiag(0)) > 1 and (cr o ff(0)) < 1, or equivalently cos 6 W < 2a < 1. In this region we have 
u) = 0 while /3 7 ^ 0, leading to a modified W ± propagator, and a free photon and a massive Z boson. 

• The remaining parts of parameter space, where (adiag(O)) > 1 and (cr Q ff(0)) > 1, or 0 < 2a < cos$w In this 
regime we have both /3 yf 0 and w / 0, which modifies the W , Z and photon propagators. Furthermore this 
region will fall apart in two separate regions III and IV due to different behaviour of the propagators (see 
Figure 6). 

5.2.1 The off-diagonal gauge bosons 

Let us first look at the behaviour of the off-diagonal bosons under the influence of the Gribov horizon. 
The propagator (112a) only contains the /3 Gribov parameter, meaning to does not need be considered here. 

As found in the previous section, this (3 is not necessary in the regime a > 1/2, due to the ghost form 
factor (odiag(O)) always being smaller than one. In this case, the off-diagonal boson propagator is simply of the 
massive type. 

In the case that a < 1/2, the relevant ghost form factor is not automatically smaller than one anymore, 
and the Gribov parameter /3 becomes necessary. The value of /3 is given by the gap equations (115), which has 
exactly the same form as in the case without electromagnetic sector. Therefore the results will also be analogous. 
As the analysis is quite involved, we just quote the results here. 

For 1/e < a < 1/2 the off-diagonal boson field has two real massive poles in its two-point function. One 
of these has a negative residue, however. This means we find one physical massive excitation, and one unphysical 
mode in this regime. When a < 1/e the two poles acquire a non-zero imaginary part and there are no poles with 
real mass-squared left. In this region the off-diagonal boson propagator is of Gribov type, and the W boson is 
completely removed from the spectrum. More details can be found in [17]. 

5.2.2 The diagonal SU(2) boson and the photon 

The two other gauge bosons — the diagonal SU(2) boson and the photon, Z M and the A M — have their 
propagators given by (112b), (112c) and (112d). Here, oj is the only of the two Gribov parameters present. 

In the regime a' > 1/2, oj is not necessary to restrict the region of integration to fi. Due to this, the 
propagators are unmodified in comparison to the perturbative case. 

In the region a 1 < 1/2 the Gribov parameter oj does become necessary, and it has to be computed by 
solving its gap equation. Due to its complexity it seems impossible to compute analytically. Therefore we turn 
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Figure 6: Left is a plot of the region a' < 1/2 (the region a' > 1/2 covers all points with higher u). In red are points where 
the polynomial in the denominator of (111b) - (Hid) has three real roots, and in blue are the points where it 
has one real and two complex conjugate roots. At the right is a slice of the phase diagram for g = 10. The region 
a > 1/2 and a ' > 1/2 is labelled I, the region a < 1/2 and a' > 1/2 is II, and the region a < 1/2 and a' < 1/2 is 
split into the regions III (polynomial in the denominator of (111b) - (llld) has three real roots, red dots in the 
diagram at the left) and IV (one real and two complex conjugate roots, blue dots in the diagram at the left). 
The dashed line separates the different regimes for off-diagonal SU(2) bosons (two real massive poles above the 
line, two complex conjugate poles below). 


to numerical methods. Once the parameter to has been (numerically) determined, we look at the propagators to 
investigate the nature of the spectrum. 

As the model under consideration depends on three dimensionless parameters ( g , g' and i///t), it is not 
possible to plot the parameter dependence of these masses in a visually comprehensible way. Therefore we limit 
ourselves to discussing the behaviour we observed. 

In region III, when there are three real poles in the full two-point function, it turns out that only the two 
of the three roots we identified have a positive residue and can correspond to physical states, being the one with 
highest and the one with lowest mass squared. The third one, the root of intermediate value, has negative residue 
and thus belongs to some negative-norm state, which cannot be physical. All three states have non-zero mass for 
non-zero values of the electromagnetic coupling g', with the lightest of the states becoming massless in the limit 
g' — > 0. In this limit we recover the behaviour found in this regime in the pure SU{2) case [17] (the Z-boson 
field having one physical and one negative-norm pole in the propagator) with a massless boson decoupled from the 
non-Abelian sector. 

In region IV there is only one state with real mass squared — the other two having complex mass squared, 
conjugate to each other — and from the partial fraction decomposition follows that it has positive residue. This 
means that, in this region, the diagonal-plus-photon sector contains one physical massive state (becoming massless 
in the limit g' —» 0), and two states that can, at best, be interpreted as confined. 

6 Conclusion 

In this work we presented results achieved during the study of Yang-Mills models, in the Landau gauge, 
coupled to a Higgs field, taking into account non-perturbative effects. More specifically, the SU (2) and SU( 2) x t/( 1) 
models were analysed in 3- and 4-dimensional Euclidean space-time, while the Higgs field was considered in its 
fundamental and adjoint representation. The non-perturbative effects were taken into account by considering the 
existence of Gribov ambiguities, or Gribov copies, in the (Landau) gauge fixing process. In order to get rid of those 
ambiguities we followed the procedure developed by Gribov in his seminal work [10], which consists in restricting 
the configuration space of the gauge field into the first Gribov region fh As found by Gribov, that restriction of the 
path integral domain leads to a modification of the gluon propagator, in a way that it is not possible to attach any 
physical particle interpretation to that quantity. The gauge field propagator develops, after the Gribov restriction, 
two complex conjugate poles, which spoils the physical interpretation of a particle. This may be interpreted as 
a sign of confinement of the gluon field. In the same sense we could observe similar modifications of the gluon 
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propagator in the probed models. In general, the poles of the gauge field propagator are functions of the parameters 
in each Yang-Mills model (including the Gribov parameter 7 and the Higgs self mass parameter v), so that we 
could identify regions in the parameter space where the gauge propagator has contributions coming from Yukawa 
type modes (with real poles) and/or contributions from Gribov like modes (with cc poles). Regions where only 
Yukawa contribution exist are called Higgs regime or Higgs “phase”. At the same time, regions where there is only 
Gribov type contributions are named confined regime or confined “phase”. Note that contributions with negative 
residue, despite being of the Yukawa type, have no physical particle interpretation as well. 

In general we could find that the Higgs regime corresponds to the region of weak coupling, i.e. small g 
and sufficiently large 17 reached in the UV regime. For that region of the parameter space (coupling constant and 
Higgs vacuum expectation value), where we expect perturbation theory to work, we do recover the perturbative 
Yang-Mills-Higgs propagators. This is an important observation, since it means that the Gribov ambiguity does 
not spoils the physical vector boson interpretation of the gauge sector, where it is relevant. On the other side, 
the confined phase corresponds to the strong coupling region in the parameter space, characterized by large values 
of g and sufficiently small v, but still keeping logarithmic divergences under control. For higher values of the 
non-Abelian coupling constant the Gribov horizon lets its influence be felt and the propagators become modified. 
In general, for the SU( 2) and for the SU( 2) x U[l) cases, we could find an intermediate region where contributions 
from physical modes (with real poles) mix with contributions from non-physical modes (with cc poles or negative 
residues). 

For the fundamental representation of the Higgs field, either in the SU( 2) or in the SU( 2) x 17(1) model, 
we could find that these two “phases” (confinement/Higgs) may be continuously connected. This feature seems to 
be in qualitative agreement with the work by Fradkin & Shenker [4]. According to [4, 7, 35] one observes a region 
in the (v,g) plane, called analyticity region, in which the Higgs and confining “phases” are connected by paths 
along which the expectation value of any local correlation function varies analytically, implying the absence of any 
discontinuity in the thermodynamical quantities. According to [I], the spectrum of the theory evolves continuously 
from one regime to the other. It turns out that these two regions are smoothly connected. In our work no rigorous 
analysis were carried out concerning the existence of any analyticity region, given the obvious complexity of the 
matter. 

The word “phase” was carefully used in this work, since the existence of a gauge invariant local operator, 
whose vev is sensitive to a phase transition, was not probed in this work. At finite temperature the use of Polyakov 
loop becomes a useful order parameter, and works on this subject point to second order phase transition for the 
SU(2) gauge theory (see [36, 37] and a more recent work [38]). 

In the adjoint representation of the Higgs field things look quite different. Besides the confining phase, in 
which the gluon propagator is of the Gribov type, our results indicate the existence of what can be called a U( 1) 
confining phase for finite values of the Higgs condensate. This is a phase in which the third component A ^ of the 
gauge field displays a propagator of the Gribov type, while the remaining off-diagonal components A“, a = 1,2, 
exhibit a propagator of the Yukawa type. Interestingly, this phase has been already detected in the lattice studies 
of the three dimensional Georgi-Glashow model [5, 6 ]. A second result is the absence of the Coulomb phase for 
finite Higgs condensate. For an infinite value of the latter, we were able to clearly reveal the existence of a massless 
photon, in agreement with the lattice suggestion, as e.g. [33]. 

That non-perturbative analysis of Yang-Mills theories coupled to a Higgs field can be extended in many 
interesting ways. For instance, a similar analysis could be done considering the all-order approach to the Gribov 
problem, known as Gribov-Zwanziger approach [39, 40, 41] (see also [11, 12] for detailed study), and also the refined 
version of it, when non-perturbative effects concerning the existence of dimension two condensates are taken into 
account [42, 43, 44, 45, 46]. Equally interesting is the analysis of these Yang-Mills models at finite temperature, 
which is already being carried out by some of the authors [38]. 
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